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Abstract 

Utilizing a QCD Coulomb gauge Hamiltonian with linear confinement spec- 
ified by lattice, we report a relativistic many-body calculation for the light 
exotic and charmed hybrid mesons. The Hamiltonian successfully describes 
both quark and gluon sectors, with vacuum and quasiparticle properties gen- 
erated by a BCS transformation and more elaborate TDA and RPA diag- 
onalizations for the meson {qq) and glueball {gg) masses. Hybrids entail a 
computationally intense relativistic three quasiparticle {qqg) calculation with 
the 9 dimensional Hamiltonian matrix elements evaluated variationally by 
Monte Carlo techniques. Our new TDA (RPA) spectrum for the nonexotic 
1 charmed (cc and ccg) system provides an explanation for the overpopu- 
lation of the observed J states. For the important 1 ^ light exotic channel 
we obtain hybrid masses above 2 GeV, in broad agreement with lattice and 
flux tube models, indicating that the recently observed resonances at 1.4 and 

1.6 GeV are of different, perhaps four quark, structure. 
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Exotic hybrids, hadrons with quantum numbers not possible in simple qq or qqq quark 
models, have been an elusive, yet signature prediction of Quantum Chromodynamics (QCD). 
It was therefore quite natural that the recent observation by the E852 collaboration |1| 
of two exotic J^^ = 1"+ states with masses 1.4 and 1.6 GeV would attract widespread 
interest. Since these states have isospin / = 1 they can not be glueballs (oddballs) and 
would initially appear to be viable hybrid meson candidates, especially since vintage bag 
model calculations predict exotic excitations with explicit gluonic degrees of freedom in 
this mass range. However, the detailed structure of these states remains uncertain since most 
contemporary theoretical studies, such as lattice gauge flux tube QCD sum rule 

0, perturbative non-relativistic QCD [|1^ and constituent models |Tl| have focused on heavy 



quark hybrids. The only two modern light quark hybrid calculations, a lattice gauge |Q and 
flux tube [§], have further compounded this uncertainty by predicting the lightest exotic 
hybrid mass to be about 2.0 GeV- significantly above the observed l""*" states. Because the 
bag model results are rather dated and lattice calculations are less accurate for light quarks 
due to extrapolation, it is important to have an additional, alternative hybrid prediction. 
The purpose of this Letter is to determine if these exotic states can indeed be interpreted as 
hybrids within a relativisitc many-body constituent approach that has successfully described 
both conventional meson [l^Jl^ and glueball systems. 



Our starting point is the QCD Coulomb gauge Hamiltonian (see for example ref. |]T5[ ) 
which we simplify to a form amenable for many body calculations 

H = j dx^^(x)(-2a-V + /3m)^(x) +Tr j d^{W -TV + 

-\ |rfxrfyp'^(x)r(|x-y|K(y). (1) 

Here and A, are the respective quark and gluon fields, = V x A'', and p°' = 
vj/tj^a^ _l_ jahcph . "Qc -g quark plus gluon color density. The current quark mass, m, 
is assigned the values, m„ = = 5 MeV and rric = 1200 MeV for the u, d and c flavors, 
respectively. Confinement and leading canonical interactions are represented by the instan- 
taneous potential, V = — — + err, with as = .2, and, a = 0.135 GeV^, as determined by the 



string tension from lattice and Regge fits. We also use a cut-off parameter A = 4 — 5 GeV to 
regularize the logarithmic divergent term in the mass gap equation. The model parameters 
a, as and A are commensurate with our previous pure quark |T2|JT3[] and gluon |T^ applica- 
tions which produced reasonable hadronic descriptions including the Regge trajectory slopes 
for the mesons (e.g. p tower) and glueballs (pomeron) [ |16| |. 

Next we proceed to the many-body diagonalizations but first perform a canonical trans- 
formation (BCS rotation) to a new quasiparticle basis 



(2) 



a^ik) = cosh Ok <(k) + sinh 9^ af{-k) 

Bcx(k) = cos y&cA(k) - A sin ^dl^{-k) 

6 6 
Dcx{-k) = cos -^4A(-k) + A sin Y^lxi^) > 



where 9fc, 9k/2 are the BCS angles, further specified below, and a{a), b{B) and d{D) are bare 
(dressed) gluon, quark and antiquark Fock operators, respectively. The indices a = 1,2. ..8 
and c = 1, 2, 3 denote color while A represents spin projection. The new field expansions are 

dk 1 



4"(x) 

nr(x) = - 

for the gluon fields and 



(27r)3 V2ZJ^ 



K(k)+af(-k)]e' 



cA 



^;fK(k)-a.f(-k)|e- 



U,xik)B,,{k) + V,x{-k)Dl{-k) 



(3) 
(4) 
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for the fermion field. The rotated Dirac spinors are given in terms of the Pauli spinors, 



Ca(k) = ^ 



VI + sin 0fc Xcx 



Vr^^sm0l a ■ kxcx 



VI - sin0fc(T ■ kxcX 



y/1 + Sin0fc XcX 



(6) 



We then find an improved, nontrivial vacuum, 1^2), by minimizing the ground state 
expectation value of the Hamiltonian variationally with respect to the BCS angles. Actually, 
the specific variational parameters are the quark gap angle, 0^, related to the BCS angle by 
tan{(f)k — dk) = Tn/k, and the gluon self-energy, Uk, satisfying Uk = ke''^®''. This generates a 



quark and gluon gap equation (equivalent to the Schwinger-Dyson equation) yielding mass 
gaps of about 100 MeV for the u/d quarks and 800 MeV for the gluon. The BCS vacuum 
contains quark and gluon condensates (Cooper pairs) in reasonable agreement with QCD 
sum rules. For more complete details consult refs. |T2|-p!^. 

We find the two body sector is adequately described by the Tamm-Dancoff approximation 
(TDA) in which a glueball is represented by Fock states g'^g^Q) and mesons by g^g^|fi). The 
notable exception is the light pseudoscalar sector, where a more sophisticated, collective 
approximation (Random Phase approximation or RPA) is needed to correctly reproduce the 
Goldstone boson nature of the pion due to spontaneous chiral symmetry breaking by our 



BCS vacuum. This is also further documented in refs. [|T2|,|T3 . 

Finally we formulate the hybrid meson as \hybrid) = q'^q^g^Q) = [[q^ 'SlQ^s'S) d'^hl^) y 
now involving quark color octet states. The resulting TDA equation for the hybrid mass M 



IS 



{hybrid\[H,q^fg^\Q) = M {hybrid\q^f g^Q) . (7) 

This projects the hybrid meson wave equation, which is pictorially represented in Fig. |l|, 
onto the three body Fock basis. Unlike our pion application, the qq pair is now in a color 
octet and the TDA is sufficient since the Hamiltonian, as well as exact QCD, does not 
conserve the chiral color octet current (we specifically calculated no difference in the more 
elaborate, chiral symmetry preserving hybrid RPA calculation, see below). The relevant 
angular momenta (am) are the q, q spins coupled to an intermediate S, and the gluon spin 
with its orbital am L_|_ (with respect to the qq cm) coupled to intermediate I. Coupling / 
with L_ (the orbital qq am) yields L which combines with 5* giving the total am J. The 
complete wavefunction in the hybrid cm has form 

i^AXA,(q+,q-) = E F/,XLsiM, |q-l) >r,"(q+)>r_-(q-) (8) 

{-ly^iL+m+l - Xg\lmi){L_m_lmi\LmL){^Xg^ - \g\Sms){-l)^~^'^{LmLSms\Jmj) , 
where q_ and q+ are the respective relative momentum of the qq pair and gluon (with 



respect to the pair cm). We then impose the transversahty condition, k- a(k) = 0, from the 
Coulomb gauge constraint which ehminates states with L+ = 1 and / = 0. For pure S waves 
the hghtest hybrid states will then have J^'" = 1+^, and 2++. These are nonexotic 

states which will mix with conventional mesons and hinder hybrid identification. For exotic 
states one P wave is necessary and we calculate the lightest corresponds to = 1 since 
the L_ = 1 excitation is energetically more expensive due to quark repulsion in the octet 
channel. This generates the exotic states 1^^, 3^^ and . 

Instead of solving the formidable TDA nonlocal equations (effectively a 12- dimensional 
problem in momentum space), we evaluate the hybrid mass variationally using an exponen- 
tial radial wavefunction for each of the two independent momentum variables. In the center 
of momentum frame the matrix elements reduce to 9-dimensional integrals that we evaluate 
numerically using the Monte Carlo code VEGAS. We then perform searches for minima on 
the energy surface in the different angular momentum channels. Our final results and key 
findings of this Letter are displayed in Figs. |^ and |[ 

Note from Fig. ^ the clear agreement between our predictions and the lattice and fiux 
tube results for both light and charmed hybrid states. This agreement sharply contrast 
with the BNL measurements which strongly suggests that the observed exotic states are 
not hybrids. To confirm our result is not an artifact of the variational method, we have 
reproduced our conventional meson and glueball exact TDA spectra to within a few percent. 
We also varied the least constrained model parameter ag, from 0.2 to 0.4. Related, we 
even performed a more extensive RPA variational calculation and took the chiral limit 
{rriu = md 0), finding only miniscule change in the hybrid mass, consistent with the 
nonconservation of the chiral color octet charge discussed above. The culmination of our 
model sensitivity study produced at most a 10 % hybrid mass variation indicated by the 
box in Fig. |^. 

Since four quark states qqqq can also have exotic quantum numbers, one can make simple 
estimates yielding exotic masses between 1 and 2 GeV for quarks in color singlet configu- 
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rations. This is consistent with a recent unitary quark model calculation which also 
concluded that the observed states are indeed predominately meson-meson resonances. 

We also calculated the lightest non-exotic hybrid (ground state) to have mass slightly 
above 2 GeV. It is interesting to note that this state has J^*" = l"*"" in contrast to certain 
heavy hybrid, quenched lattice ground state results which find near degeneracies among 
several negative parity states. This was first noted by ref. |]l8l who used a similar constituent 



model in the heavy (static) quark limit. However, the heavy quark lattice simulations do 
not include quark spin and until this is included the ground state hybrid quantum numbers, 



as well as related level ordering, remains unclear [|19|. Further, another lattice calculation 
studying orbital hybrid excitations concluded that the ground state quantum numbers were 
likely to be l"*"" in agreement with our work. This study noted that the degeneracies in the 
static framework would be broken by quark spin-orbit effects which shift the exotic levels 
and also mix nonexotic and conventional qq states. They were also able to infer the level 
splitting, yielding the above conclusion, using propogating quarks on the lattice. There 
is less, but still some, uncertainty in level order within our model as well since our spin 
dependent interaction is incomplete. Although further level ordering study is necessary, the 
thrust of this work is the exotic hybrid and improved spin effects will not alter our conclusion 
that the lightest is above 2 GeV. 

In Fig. we compare our full model spectrum to data for the well studied, believed to be 
gluon rich, 1 , J/ip system to provide an explanation for the anomalous overpopulation of 



observed states ||2^ with respect to quark model predictions. Whereas previous constituent 



calculations, using only S waves, could only account for 3 of the known 6 charmonium 



levels, we now predict 7 cc states in addition to 4 ccg hybrids. Further, ref. pO[ lists an 
additional charmonium level iIj{383Q) assigned J^*" = 2 which also agrees well with our D 
wave prediction (not shown). Notice that by simply including D waves we have resolved the 
"overpopulation" problem. In general all of the states, both cc and ccg, will mix and 
a more elaborate calculation is in progress. However, our current result is already sufficient 
to conclude that simple level counting (density of states) arguments will probably not be 

6 



effective in identifying charmed hybrid states. 

Finally, we mention a novel color octet effect leading to an isospin splitting since it only 
affects the / = states. This is the annihilation process depicted in Fig. |I] corresponding to 

^ 9 ~^ for the L_ = 0, spin aligned color octet quark pair. Octet quarkonium is the 
QCD analogue to ortho positronium and the annihilation interaction raises all / = light 
hybrid states by roughly 300 MeV when the q — q spins are aligned {S = 1). 

Summarizing, our large-scale diagonalizations of an effective Coulomb gauge Hamiltonian 
provide a reasonable, comprehensive description of the meson, lattice glueball and lattice 
hybrid meson spectra. Further, our composite J/tp spectrum is now also in much better 
agreement with data, especially in terms of density of states. It is important to note that 
our quark/gluon unified approach essentially entails only one pre-determined dynamical 
parameter. Finally, and perhaps most significant, our reaffirmation of lattice and flux tube 
1 ^ masses indicates that the recently observed exotic states are not hybrids. Based upon 
preliminary estimates and other independent studies it is more likely that these resonances 
are four quark states and more rigorous, higher quark Fock state calculations are in progress. 

We thank NERSC for providing Cray J-90 CPU time. F. L. E. acknowledges SURA- 
Jefferson Lab for a graduate fellowship. This work was partially supported by grants DOE 
DE-FG02-97ER41048 and NSF INT-9807009. 
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FIGURES 
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FIG. 1. One (dot) and two-body (waves) TDA matrix element for hybrid mesons with one 
constituent gluon. Note the qq annihilation two-body diagram. 



10 



Lattice 



TDA 



Flux Tube Bag Model Experiment 



5000 



4000 




Mass 
(MeV) 



3000 



Charm Sector 
Light Sector 



PC - + 
J = 1 



2000 



1000 



FIG. 2. Comparison of exotic 1 ^ u/d and c hybrid masses with alternative theories and data. 
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FIG. 3. TDA theory for conventional (cc) and hybrid {ccg) states compared to the observed 
1— J/ip spectrum (PDG2000) from ref. [20] 
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